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Abstract. We obtain a purely local characterisation that singles out the Majumdar- 
Papapetrou class, the near-horizon Bertotti-Robinson geometry and the Reissner- 
Nordstrom exterior solution, together with its plane and hyperbolic counterparts, 
among the static electrovacuum spacetimes. These five classes are found to form 
the whole set of static Einstein-Maxwell fields without sources and conformally flat 
space of orbits, this is, the conformastat electrovacuum spacetimes. The main 
part of the proof consists in showing that a functional relationship between the 
gravitational and electromagnetic potentials must always exist. The classification 
procedure provides also an improved characterisation of Majumdar-Papapetrou, by 
only requiring a conformally flat space of orbits with a vanishing Ricci scalar of 
the usual conveniently rescaled 3-metric. A simple global consideration allows us to 
state that the asymptotically flat subset of the Majumdar-Papapetrou class and the 
Reissner-Nordstrom exterior solution are the only asymptotically flat conformastat 
electrovacuum spacetimes. 
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1. Introduction 



The (standard) Majumdar-Papapetrou and Reissner-Nordstrom metrics are known to 
describe, under rather general conditions, the exterior geometries of the static charged 
black holes, as shown by the recent uniqueness theorems (see [1] and references therein). 
The aim of this work is to provide an essentially local characterisation for the Majumdar- 
Papapetrou class and the Reissner-Nordstrom exterior solutions. Local characterisations 
are important, not only for being essential ingredients for the improvement of the global 
charaterisations of black holes provided by the uniqueness theorems, but also for a 
better understanding of the solutions and its potential use in stability problems. We 
first find a purely local uniqueness result that characterises Majumdar-Papapetrou, the 
near-horizon geometry and Reissner-Nordstrom, together with its plane and hyperbolic 
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counterparts, among the (strict) static electrovacuum spacetimes. Global considerations 
can then be used to restrict the set conveniently. 

The known characterisations of the different families of black hole metrics vary 
from the (purely) local, the "essentially" local and those of global nature in quite a 
gradual manner. Whether or not some global property is preferable to some stronger 
local constraints is not clear 'a priori' (see e.g. the discussion about some possible Kerr 
characterisations included in |2]). Nevertheless, it is always convenient to try to minimise 
the number of constraints present in any given local characterisation. Take the different 
characterisations we have for the Schwarzschild metric. Although Birkhoff's theorem 
constitutes a nice and purely local characterisation, it seems of no use in the uniqueness 
theorems. Another purely local characterization which involves only the Weyl tensor and 
the metric itself is given in [3]. A more convenient characterisation ingredient appears 
to be the conformal flatness of the hypersurfaces of constant static time (conformastat), 
since that constitutes a crucial step in the uniqueness theorems as they stand now. 
This is, in fact, not purely local, since this characterises Schwarzschild among the static 
and asymptotically flat vacuum spacetimes. Indeed, conformastat vacuum spacetimes 
comprise three |1] (see also [5]) out of the seven families that constitute the whole set of 
degenerate (type D) static vacuum spacetimes [6l[7]. These three families correspond to 
the Schwarzschild solution together with its plane and hyperbolic counterparts (Class 
A in Table 18.2 in [7]). The Schwarzschild solution can be singled out by requiring 
asymptotic flatness. This is a simple global consideration, and in this respect one may 
think of this as being an essentially local characterisation. 

A natural step to follow is the generalisation of the above to static charged black 
holes. Indeed, global arguments in the uniqueness theorems establish, again, conformal 
flatness of the hypersurfaces of constant static time [Sj. It is important to note, however, 
that the same global arguments also imply that the gradients of the gravitational and 
the electromagnetic potentials are aligned, or in other words, that the potentials are 
functionally related. It is these two facts, together with a "non-degeneracy" restriction 
and asymptotic flatness, that lead eventually to spherical symmetry and thus to the 
standard uniqueness results for the non-extreme Reissner-Nordstrom black hole (see 
e.g. [8]). One question we address in this paper is up to which extent the alignment and 
the "non-degeneracy" properties can be relaxed in a local characterisation. We believe 
this may be of use on the improvement of the recent uniqueness theorems of (multi) 
black holes (see [1] and references), since we provide an essentially local uniqueness result 
for the static charged black hole solutions, binding together the Majumdar-Papapetrou 
to the exterior Reissner-Nordstrom solution. 

The characterisation we present here come by finding the complete solution 
of the Einstein-Maxwell field equations without sources for static spacetimes with 
a conformally flat space of orbits. The only extra assumption made is that the 
electromagnetic field inherits the symmetry, so that it is also stationary. We call such 
solutions conformastat electrovacuum spacetimes. Note that we take the definition in 
[7] as standard, following the original terminology by Synge [9]: conformastationary are 
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those stationary spacetimes with a conformally flat space of orbits and the conformastat 
comprise the static subset. 

Conformal flatness corresponds to the vanishing of the Cotton tensor associated 
to the induced metric on the space of orbits. A general study of conformastationary 
spacetimes would follow then an analogous path to the chatacterisation of the Kerr 
and Kerr-Newman families of black holes among the stationary solutions. In the Kerr 
case the crucial local property is the vanishing of the complex Simon tensor [10] , which 
generalises the Cotton tensor on the space of orbits. The characterisation of the Kerr 
metric in [10] comes as a result of the equivalence of the multipole structure of Kerr with 
that of an asymptotically flat end with vanishing Simon tensor. The first objection to 
this characterisation is, precisely, that the isometry with Kerr is only established in some 
neighbourhood of infinity, and hence the extension of this isometry to the whole (strict) 
stationary region cannot be ensured yet. This motivates, in fact, the search for improved 
local characterisations, since the problem of the extension of the isometries to whole 
(strict) stationary regions may be fixed by exploiting the local characterisations to their 
full extent. Indeed Perjes found [13] that the most general metric with vanishing Simon 
tensor depends only on a few parameters, and thus showed that the asymptotically 
flatness condition in the characterization of Kerr is only necessary in order to fix the value 
of some constants. In this paper we thus follow an analogous aim, since we exhaust the 
implications of the vanishing of the Cotton tensor in the static electrovacuum problem. 

The second drawback the characterisation of Kerr in [10] faces is that, by 
construction, it is not valid within the ergosphere. To address this problem Mars [Hj |2] 
managed to improve that characterisation and include the ergosphere by constructing 
the so called Mars-Simon tensor, this time relative to the spacetime. The Kerr 
characterisation in [11] (see also Theorem 1 in [2]) is essentially local, since the vanishing 
of the Mars-Simon tensor produces a family of vacuum solutions depending on two 
complex constants, only to be fixed by some simple global consideration. On the other 
hand, in [2] , Mars provided a characterisation with a much weaker local condition, using 
more effectively the asymptotic flatness. The work in [11] has been extended recently by 
Wong in [12] , by providing a couple of extended characterisations for the Kerr-Newman 
family, the first being purely local. 

The main assumptions inherent to the spacetime characterisations of the Kerr- 
Newman family have two crucial direct implications. The first is the degeneration of the 
Weyl tensor (type D), and the second is the existence of a functional relationship of the 
gravitational and electromagnetic potentials in the static case. None of these restrictions 
are taken as assumptions in the present work. Not imposing any restriction on the 
Petrov type is important in the static otherwise the Majumdar-Papapetrou 

class would not be taken into cosideration. On the other hand, the key result in the 
present paper that leads to the complete solution of the conformastat electrovacuum 
problem is precisely that the aligment of the gradients of the potentials is necessary. In 
this sense, in the static case the results found here generalise completely those in [T2] . 
Furthermore, these results suggest that the known local Kerr-Newman characterisations 
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may be improved by relaxing some of the requirements involved. 

The vanishing of the Cotton tensor in the stationary vacuum problem was dealt 
with in a series of three papers by Lukasz et al. in [11] and Perjes in [151 IB] (see 
also [H]). They found the whole set of conformastationarj||| vacuum spacetimes. In a 
first paper [H] the solutions possessing a functional relationship between the real and 
imaginary parts of the Ernst potential S were found to consist of three bi-parametric 
families of solutions generated from the three conformastat vacuum solutions (Class A) 
by the Ehlers transformation. In [15] , using the purposely defined "Ernst coordinates" , 
Perjes found that solutions with functionally independent real and imaginary parts of S 
necessarily admit a spacelike isometr}j§|, to conclude in [16] (see also [H]) that this set 
of solutions is empty. Therefore, all conformastationary vacuum spacetimes belong to 
the three families presented in [14], which can be thought as the NUT- type extensions 
of Schwarzschild and its plane and hyperbolic counterparts. 

The plan of this paper is analogous. We start in Sections |2] and |3] by showing how 
the conformastat electrocavuum problem and the conformastationary vacuum problem 
can be treated within a common framework by using a suitable notation. The motivation 
is to use the previous works [151 IlHl IE] as a guide, and additionally, to recover 
those results. In Section H] we prove the key result: the conformastat electrovacuum 
spacetimes necessarily contain a functional relationship between the gravitational and 
electromagnetic potentials. Regarding the use of the procedures in references [151 [16] two 
points must be stressed. Firstly, the "common" proof needs at many stages a different 
approach, since the variables involved in the general case are not necessarily complex, 
and thus the positivity of some products cannot be used. Secondly, the final stages in 
the proof differ from those in |16] and fix some errors found in [17]. Anyway, to ease 
the comparison with these works we have kept the same notation whenever possible. 

In the second part. Section [5] is devoted to complete the study of conformastat 
electrovacuum spacetimes by classifying and exploiting the necessary functional 
relationship between the gravitational and electromagnetic potentials. We find that 
all conformastat electrovacuum spacetimes either belong to the Majumdar-Papapetrou 
class or correspond to either the Bertotti- Robins on solution or the exterior Reissner- 
Nordstrom solution toghether with its plane and hyperbolic counterparts. Furthermore, 
the procedure used for the classification provides an improved characterisation of the 
Majumdar-Papapetrou class. This is known to be the class of static electrovacuum 
spacetimes such that the usual rescaled induced metric in the space of orbits is flat. 
Here we find that one only needs to ask that metric to be conformally fiat and with 
vanishing Ricci scalar. 

The main result constitutes then a completely local characterisation of the 
static and charged (multi) black hole solutions, plus the "non-standard" Majumdar- 
Papapetrou solutions, the near-horizon geometry (Bertotti-Robinson) and the plane 

I Let us note that they refer to conformastationary spacetimes simply as "conformastat" . 

§ The authors talk of an "axial" symmetry, but no global property is involved in the result at this 

point. 
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and hyperbolic counterparts of the exterior Reissner-Nordstrom. A simple global 
consideration can be used now to single out the black hole solutions. The essentially 
local characterisation is thus that the conformastat electrovacuum asymptotically flat 
spacetimes are isometric either to the asymptotically flat subset of the Majumdar- 
Papapetrou class or the Reissner-Nordstrom static exterior. 

Sign conventions differ from those in [TU [151 CSj and follow those in [7] : the metric 
has signature (—,+,+,+) and the Riemann tensor is defined so that 2'W[qVp]W^ = 
R^jaisW^- Greek indices refer to the spacetime and Latin indices to the three-dimensional 
space of orbits. Units are chosen so that G = c = 1. 

2. Conformastationary spacetimes 

A stationary spacetime {A4,gnu) is locally defined by the existence of a timelike 
Killing vector field whose space of orbits invariantly determines a differentiable 3- 
dimensional Riemannian manifold S3. Local coordinates exist for which ^'^ = dt 

and such that the line-element can be cast as [7j 

ds^ = -e^^{dt + AJx^Y + e~^^habdx^dx\ (1) 

where U, A^ and h^h do not depend on t. Applying the usual projection formalism 
[m [7j we will think of \J function on S3, Aa l-form belonging to T*S3 and 
hah as a metric on S3. Once these three objects are given, the local geometry of the 
stationary spacetime {M.,g^v) is fully specified by using ([T]). Let us, from now on, 
endow S3 with the metric h^h and use the first latin indices a,b, . . . for objects defined 
on (S3, /lab). A conformastationary spacetime is a stationary spacetime whose space of 
orbits (S3, hab) is conformally flat [7]. Thence, in a conformastationary spacetime there 
exist coordinates {x, y, z} in which habdx°'dx'' = e^^^''^''^\dx^ + dy"^ + dz^). The intrinsic 
characterisation of a conformally flat 3-space is the vanishing of the Cotton tensor Cabc 
[H [19], or equivalently, the York tensor density [20], defined as Yg^ = fj^'^^Cabc, where 
fjabc denotes the volume form of (S3,/iaf,), which satisfies Yae = Yea and Ya" = 0. More 
expliclitly, S3 is conformally flat if and only if 

Ya' = fj'^' (2VMba - \habVcRj = 0, (2) 

where Rab and V denote the Ricci tensor and covariant derivative relative to hab- 

Conformastat spacetimes are those conformastationary spacetimes which are, in 
fact, static. In this context, a static spacetime is thus characterised by Aa = 0. 

2.1. Electrovacuum field equations 

Let us first fix one basic assumption and some notation. First, we will restrict ourselves 
to Maxwell fields Fq,^ in {M,gfj,u) which inherit the stationary symmetry, i.e. for which 
C^F = 0. The Einstein- Maxwell equations outside the sources imply (locally, at least) 
the existence of two complex scalars, $(a:°) the electromagnetic potential, and S{x'^) the 
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Ernst potential. These two potentials in (S3, hat) satisfy the so-called Ernst-Maxwell 
equations, 

V''Ha + ^G-H-^G-H = 0, (3) 

V^Ga -H-H-{G-G)-G = 0, (4) 

where Ha = {Re£ + and Ga = l/2{Re£ + $$)-H^,a + 2$<l>,a) and the dot 

denotes the scalar product. It will be convenient for later to note two identities that Ha 
and Ga satisfy: dH = H A ReG and dG = GAG + HAH. These relations are, in fact, 
the integrability conditions for the two potentials. 

The rest of the Einstein- Maxwell equations without sources reduce to the following 
problem for hat 

Rab = GflGb + GaGf, — {HaHb + HaHb). (5) 

Once hab is known, the geometry and electromagnetic field are recovered from the 
complex potentials. The metric function U and the 1-form Aa are determined by the 
relations 

e^^ = Re^ + $$, {dA)ab = 2e-^'' f^abJ^mG'' , 

taking into account that the freedom in the determination of Aa corresponds to a 
transformation of the time coordinate of the form t — )■ [7]. The electromagnetic 

field, conveniently described by the self dual 2-form J^^^ 

= F + i* F, 

where * stands for the Hodge dual in {M.,g^y), i.e. *Fai3 = \F^^rjfj,ua.i3, is thus recovered 
by 

:F = -e-^ [H Ai + i*{H Ai)], 

where the 1-form H^ in {M.,g^u) is given, in coordinates adapted to the Killing ([T]), 
by H^ = {0,Ha). Note that = — e^'^(l, The real and imaginary parts of H^ 
correspond to the electric and magnetic fields with respect to the observer defined by 
u = e~^^, this is H^ = + iB^ = F^yU'^ . For completeness, let us note that the 
intrinsic definition of G^ in (A^, g^y) is given by G^ = u^iV + i* VyU^) and its real 
and imaginary parts G^ = + correspond to the acceleration and twist vectors 
of the congruence u. 

2.2. Vacuum and electro-magnetostatic cases 

The stationary vacuum case is characterised by $ = 0, so that Ha = and hence ([5]) 
specialises to 

Rab = GaGh + GaGb 

and the Ernst-Maxwell equations reduce to 
V"Ga - (G - G) ■ G = 0. 
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The integrabihty condition that Ga satisfies reads simply dG = G A G. 

The static case is characterised by Ga — Ga{= 21mG'a) = 0. Well known fact is that 
the conditions for Ha and Ga and the field equations yield dG = (in fact Ga = U^a), 
and Ha = e~'^^^Ha for some constant 9 (see e.g. [21]). Let us now define the vector 
Xa = e~^^Ha, which is real by construction and related to the electric and magnetic 
static fields by Ea = cosOXa and Ba = sin6'Xa. Instead of working with the complex 
$ let us consider the real potential = e~*^$, so that Xa = e'^'^^a- We are thus left 
with two real vectors: Ga and Xa- 

The stationary vacuum and the static electrovacuum cases are known to have an 
analogous structure, although they are inequivalent (see e.g. Chapter 34 in [7]). The 
analogy has been used previously in the literature in a more or less implicit manner (see 
e.g. [22]). The fact that the two problems are inequivalent comes most notably from 
the signature of the potential spaces, which differ in the two cases. Despite this, one can 
make the analogy explicit, and useful in the present study, incorporating that change of 
signature by making use of a hyperbolic-complex or motor number construction, based 
on the real Clifford algebra Gii^Q{R) (see e.g. [23] and references therein), for the static 
electrovacuum problem. We call j the hyperbolic imaginary unit, which satisfies = 1, 
and denote the conjugate operation by j = —j. Note that C£o,i(-R) is isomorphic to the 
field of complex numbers, in which i is the elliptic imaginary unit. 

We are now ready to define 

J:a = l{Ga+jXa), 

in terms of which the Eintein-Maxwell and Ernst-Maxwell equations read 

Rab = ^^ah + ^a^b), (6) 
V'^S^ - (S - S) ■ S = 0, (7) 

and the identities for Ga and Xa reduce to dS = S A S. 
3. A common framework 

For the sake of completeness and to allow us to use the techniques and some results of 
previous works on conformastationary vacuum spacetimes [HI [15] , we set up a common 
and more general problem using a common notation. 

Let us denote by l any of both the complex i and the hypercomplex j, so that 
l"^ = ±1 accordingly, and the general conjugation by ~, so that T = —l stands for either 
i = —i or j = —j . Any object of the form F = f + tg will be called a composed object, 
and 9fJ(F) and '^{F) will denote its real and imaginary parts. 

Consider now a composed vector field and a real metric hab which satisfy the 



system of equations 

Rab = N{yay, + yayb), (8) 

v^ya -{y-y)-y = o, (9) 

dy = yAy. (10) 
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One could regard this problem at the level of the potentials, but for our purposes 
it suffices to set up the problem for the vectors and thence include the integrability 
conditions as equations. The vacuum case is recovered by taking N = 1, ya = Ga, 
a complex 1-form, and the conjugate being the complex conjugate. The static case 
corresponds to N = 4, y^a = S^, a j-l-form and the conjugate being the j-conjugation. 
Note that in both cases the right-hand side of the equation for Rat is, as it should, a 
real quantity, whereas the equations ([9]) and ( 1T0|) yield two real equations each. 

3.1. Conformastationarity 

Conformastationarity follows by the vanishing of the York tensor density of hat, this is, 
by applying equation ([2]) to the Ricci tensor as expressed in ([H]). Before writing down 
the explicit expressions, let us introduce a very convenient vector (see [H]) [Jj] 

where -k denotes the Hodge-dual in (T,3,hab), i-e. La = y'^y^fihca- By construction we 
have La = —La and L ■ y = L ■ y = 0. Note also that L = -kdy. Let us stress the fact 
that since La = —La, La is imaginary and thus l^{L ■ L) > 0. Introducing ([H]) into ([2]) 
one obtains the real equation 

^Ya' = {ya - yaW + ^''l5^^> VeJ^a + J^feVeJ^a) " ^ WVe(3^ " 3^) = 0. (11) 

Since Ya"' = this equation contains at most 5 independent components. We will exploit 
the consequences of those equations later. 

Two very different situations arise in the study of the system of equations composed 
by (E]), ( fTOl) . (IH]) and ( fTTI) . for ya and hab- the class of solutions for which 7^ and 
those for which La = 0. Nevertheless, before entering into the study of these two cases 
one has to consider the case y ■ y = 0. In the static case ya = is j-composed and 
S ■ S = implies, in particular, G ■ G + X ■ X = 0, which clearly leaves us only with 
the trivial case Ga = Xa = 0. However, in the vacuum case = Ga is complex and 
one can have, in principle, fields for which G ■ G = 0. The study of these null fields 
was performed in [H], where it was proven that no null coformastationary vacuum 
spacetimes exist apart from the trivial case of fiat spacetime. In the following we will 
therefore take y ■ y ^ without loss of generality. 

4. The class La ^ 

In this section we prove that the class 7^ is empty in two steps. We first show 
that if La 7^ there must be an additional isometry, and then that the existence of that 
isometry implies the non-existence of solutions with La 7^ 0. 

II Although wc have kept the notation as close as possible to that used in [T^, the vector L defined 
here differs by a multiplicative i. 
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h 



ab 



LaU + {y ■ y)yayb + {y ■ y)yayb - {y ■ y){yayb + yayb) , (12) 



Let us take the basis {La,ya,ya}- (Note that the associated basis for the real 
tangent vector space is composed by iLa, 3^a + 3^a and — 3^a)-) The metric hab 
expressed in this basis reads 
1 

L-L 

where L ■ L = {y ■ y){y ■ y) - {y ■ yy. Since ^ we have l'^{L ■ L) > 0. 

Using the obvious notation Yj^yb by Fq-^, etc., the tracefree property of Ya'^ 
translates onto 

= -{y ■ y)Yf - {y ■ y)Yyy + {y ■ + f/). 

Together with the use of the conjugate operation, this allows us to keep all the 
information contained in Yj' in only three components: Yy-^, Yy-^ and Y^-^ . (Note 
that Yy^ = Yy-^.) The corresponding three equations in f|TT]) . from where the five real 
independent equations eventually follow, read 

L'Vciy ■ 3^) = 0, (13) 
LAdL = 0, (14) 

L^L"v,3^, - ^v'^'yeLbV.iy ■y) = o. (15) 



The interpretation of the equations f|T3l) and f|T4|) is straightforward. Equation ([T] 
states that La is hypersurface orthogonal, i.e. integrable. Equation f|T3|) implies that 
the product 3^ ■ ^ is constant along La- In the static case this translates to the fact that 
the two scalars + and G ■ X are constant along the direction orthogonal to the 
planes spanned by Ga and Xa- 

4-1. The additional isometry 



In this subsection (together with Appendix A ) we prove that the above equations ( [T3 



( !T4|) and ( !T5|) . together with the Ricci equations ([8]) and the integrability condition (fTOjl 
imply the existence of a further isometry along La- 

Since La is integrable f|T^ and imaginary, there exist two real functions x(^") and 
ip{x"') such that 

L = Lxdip. (16) 

The function ip cannot be constant precisely because La 7^ 0, and we can also take x > 
without loss of generality. The integrability equations (fTOj) imply, in turn, the existence 
of two further real functions, encoded in the composed potential a{x"') so thatl^ 

y = -^da. (17) 
(7 + cr 

The main idea is to use the three potentials o", a and (f, as coordinates. In the 
vacuum (complex) case [15] these particularise to the so-called Ernst coordinates. The 

% A simple inspection shows the relationship of a with the original potentials. In the static case one 
has (T = ^{a^ + j"^): whereas in the vacuum case one recovers the usual Ernst complex potential in 
vacuum (t{= £) = e^^ + iVl. 
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independence of (p and a is ensured by the orthogonality of La and 3^^. Let us label this 
coordinate system as 

= a, x'^ = a, x^ = (f. 

The real coordinates and manifold related quantities can always be recovered by the 
obvious linear transformations to the coordinates a + a and i(cr — a). There exists a 
freedom in choosing ip, since La is invariant under the transformation 

^^/(^), X-x(l)", (18) 

for any smooth function / with non-vanishing derivative. This freedom will be only 
used in the last step of the proof (see Appendix A). 

The form of the metric hab in these coordinates follows directly from (fT2ll together 
with f lT6|) and ( !T7|) . With the help of a shorter notation for the products 

a = {a + a)\y -y), 7 = 5 = (a + a)\y ■ y), 

p = {a + d)\y-y), 

(where note that a is composed, we denote by 7 its conjugate and /3 = /3 is real) together 
with the auxiliary real functions p, which essentially substitutes and D defined by 

2 2 

D = a^-(3^ = {a + ay{L-L), = ^ = t^D^^a + a^x^ > 0, (19) 

the line-element reads 

habdx^dx^ = [-f{dx^f - 2f3dx^dx^ + a{dx^f + p^D{dx^f'\ . (20) 

Since we are dealing with 3^-3^ 7^ 0, a cannot vanish, and in the complex case one 
thus readily has that 07 > because a7 = aa. But in the j'-composed case this is not 
ensured a priori. Nevertheless, the real function D satisfies i^D > by construction, 
which in the j-composed case translates onto -D > and therefore 07 > necessarily. 
To sum up, in any case we have 

> 0. 

Let us also remark that det h = p^D"^ in these composed coordinates. On the other 
hand, given (fT6ll and (ITTj) together with the definition of L^, the volume element is 
fixed by 17123 = LpD^^y/i?D. We will take the metric to be determined by the four real 
unknown functions encoded in a, (3 and p. Without loss of generality we take p > 0. 

It only remains to write equations iQ, ([8]) plus f|T3l) and ffTSjl in this coordinate 
system. Since 



p(cr + ct)2 

equation (fT3|) holds iff a and 7 are functions of a and a only. With this information at 
hand equation (|9]) translates onto 

{ad, + (3d2) In + d,a + d^l^ = (21) 
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while (fT^ reads 

2{adi + (3d2)\np + d2/3 ^ = 0. (22) 

cr + a 

The components of the equation for the Ricci tensor ([8]) yield the four independent 
composed equations 

^ ^ ^ ^ N 

RU = Rl3 = -R33 = 0, = 7 , ^.r, , (23) 

(cr + a)^ 

which encode the six real equations, due to the fact that R22 = Ru, R23 = Ris- Note 
that R12 and R33 are real. 

The first consequence the integrability conditions of the system (pT|) -(l23l) provide 
is the following: 

Proposition 1 Any solution of the system of equations ^21\ ) . ^2^) and ( f^) for l^D > 

necessarily has /33(= 8^(3) = 0. 



To ease the reading the proof is left to Appendix A 



The only remaining function in the line-element (!20|) which may still depend on 
is p. But this cannot be the case due to 0221) . Assuming that a solution to (1221) exists, 
integration of 7 ([22]) - yields 



lnp= r /3/3i-7/32 + 2/3^^rfxM +lnpo. 
JO 219 \ + a 



The integral does not depend on x^, but the arbitrary term po depends on x^, in 
principle. However, this term can be eliminated by using remaining freedom in choosing 
the coordinates fflSl) . given by a transformation x^' = x^' {x^). 

This completes the proof of the exitence of an additional spacelike isometry 
whenever y Ay and 3^ ■ 3^ 7^ 0. 

4^.2. The class La is empty 



In Appendix B we prove the following result: 

Proposition 2 There is no solution of the system of equations f21\) . [2^ and [23^] for 
functions depending on a and a only, with l^D > 0, for N ^ 16, —2, 8/5. ■ 

Since we are interested only in the cases = 1 and = 4 we do not investigate further 
the compatibility of ([2l]), (El]) and ([23]) for the special cases A^ = 16, -2, 8/5. 

This proposition thus states that the class La ^ with and additional isometry is 
empty. Combined with Proposition [1] this finally implies that the full class 7^ is 
empty. 

We are thus only left with La = necessarily. This means that and are 
parallel, and dy = by (flO)) . Therefore is a gradient of some composed potential 
whose real and imaginary parts are functionally dependent. 

In particular, on the one hand we have thus recovered the result found in the series 
of papers [H [151 [16] (see also [IT]): 
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Theorem 1 Conformastationary vacuum spacetimes are always characterised by a 

functional relation between the potentials U and Vt. ■ 

On the other hand, in the stationary electrovacuum case we have thus proven: 

Theorem 2 Conformastat electrovacuum spacetimes are always characterised by a 

functional relation between the potentials U and $. ■ 



5. The complete solution of the conformastat electrovacuum problem 

In the conformastationary vacuum case the complete solution is thus given by those 
spacetimes for which U = U{Q). This was studied in [T3]. The solution consists of three 
explicit bi-parametric families of line-elements, as described in the Introduction. We 
refer to [H] for the explicit form of the line-elements. 

In the following we focus on the the static case. From the above results we know 
we only have to look for solutions for which 

U = f/(^). 

This is a well known ansatz used to find electro (-magneto) static solutions as described 
in [7], Section 18.6.3. Our work consists on finding all the conformastat solutions among 
this class. 

The divergence equation ([9]) firstly fixes the functional relationship to b^ (see e.g. 

0) 

e^^ = 1 - 2c^ + 

for an arbitrary constant c, which can be rewritten in parametric form in terms of an 
auxiliary function V as 

0^ = 1: = c-l/V, e^^ = (24) 

> 1 : * = c - - 1 coth V, e^^ = (c^ - 1) sinh"^ V, (25) 
c2<l: ^ = c - Vl - c2 cot V, e^^ = (1 -c2)sin-2l/, (26) 

and secondly implies 

VV = 

in all cases. The Ricci equations ([H]) reduce now to 



c' = l: 


Rab = 0, 


(27) 


c2>l: 


Rab = '^VaVb, 


(28) 


c'<l: 


Rab = —'^VaVb- 


(29) 



The remaining equation that hab and Va have to satisfy corresponds to the conformal 
fiatness of hab, and is encoded in (fTTl) . 



+ The relationship one obtains is in fact e^^ = b + a'^ + for arbitrary constants a and b. The 
constant b can be rescaled by using the freedom ^ 4* + const, (if 5 < 0) or a rescahng of the t 
coordinate (if 6 > 0). 
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Let us stress the fact that either case > 1 or < 1 constitutes a more general 
problem for hat than the problem for the conformally fiat 3-metric one encounters in the 
black hole (global) uniqueness theorems (see e.g. [8]). In the uniqueness theorems for 
charged black holes one establishes from global considerations (using the positive mass 
theorem) not only the conformal flatness of the 3-metric and that the potentials are 
functionally related, but also that the conformal factor depends only on the potential. 
Since the conformal factor is not fixed a priori in the present study, we cannot use the 
usual results found in the uniqueness theorems. Instead we follow the procedure used 
by Das in [1] in the obtaining of the static vacuum solutions. 

5.1. Case = 1 

Equation (!27|) does not involve Va and simply implies that hat must be fiat, which in turn 
renders (fTTl) to be automatically satisfied. This is the well known Majumdar-Papapetrou 
class of solutions [7]. Given any solution V of the Laplace equation V^V = in flat 
3-space, the metric of the corresponding member of the Majumdar-Papapetrou class is 
found by using ( 124|) . and thus reads 

ds^ = -^dt^ + V^idx^ + dy'^ + dz^), 

while the electromagnetic potential $ = e*^\E', after a trivial shift, is given by 

5.2. Case > 1 

In this case we are looking for solutions {hah, with \4 = 1/^ 7^ of the system 

Rab = (30) 
VV = 0, (31) 
4V[,V,]K-/ia[bV,](y-r) = 0, (32) 

where the latter stands for ([2]). 

Because of V^V = and V^a 7^ 0, local coordinates {x, y^} with A = 2, 3 can be 
chosen so that V = and also such that {y^} span the surfaces 6*2 orthogonal to Va- In 
these coordinates adapted to Va equation ( 132|) implies the following form of the metric 

habdz'^dx^ = W\x)dx^ + W{x)VtABdy^dy^ , (33) 

where W{x) is an arbitrary positive function and VLab is a Riemannian metric 
on 5*2, depending only on {y^}. The imposition of fISUl) leads to an equation for W{x) 
whose solution reads 

W = (Ae^ + 5e-^')"^ (34) 

with constants A and B. 

It only remains to see that the surfaces {S2,^ab) are of constant curvature. Let 
us consider the unit normal to 5*2, Ua = WVa, and two vectors tangent to S2, ca"', this 
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1 i^AC^BD — ^AD^Bc) 



is naCA"" = 0, such that WQab = ^A^'^B^hah- The second fundamental form of 5*2 in 
S3 thus reads Kab = ^a^^b^^ a^-b = W / {2W)VLab- On the other hand, taking into 
account the identity between the Riemann and the Ricci tensors in a 3-dimensional 
space, equation (!30!) is used to obtain the following expression of the Riemann tensor of 
kab projected on S2 

Rabcd^A^'^B^^^c'^^D'^ = ^AC^BD — ^AD^BC- 

This expression is then introduced into the Gauss equation in order to obtain the 
Riemann tensor for WVLab on 5*2, 

^"^^^ Rabcd = RabcdeA^'eB^'ec^eD'^ + KacKbd — KadKbc = 

_ fw^ ■ 

= - AABW{QaC^BD - ^AD^BC)- 

The Riemann tensor for VLab on 5*2 thus reads 

^""^RaBCD = W-'^'^^'^RabCD = -^AB{nAC^BD - ^AD^BC). 

Therefore (5*2, ^ab) is a surface of constant curvature —4AB. In principle, three different 
possibilities arise: (i) AB = 0, (ii) AB < and (iii) AB > 0. 

5.2.1. Case (i) This case is characterised by a flat VLab- Coordinates {d,^} can 
therefore be chosen such that 

By changing x — > —x if necessary, we can take 5 = and A ^ without loss of 
generality, so that W = A~'^e~'^^. The line-element and electromagnetic potential are 
now obtained by introducing this into fl55]) and using fl2^ . By performing the change 
iln(6/(r + b)) with b = {AA^y/c^ - 1)"^ > 0, together with r = 2 Vc^ - 1 1, which 



X 



induces the rescaling ^ — )■ /(2\/c'^ — 1), the line-element can be finally cast as 

ds' = -i!:±^dr2 + -^dr' + r\d^' + V), (35) 

after a further convenient rescaling of (f. The electromagnetic potential, after a trivial 
shift, reads 

r 

Note that the only restriction of the ranges of the coordinates is on r. Since we have 
taken 6 > we are left with two different ranges, — 6 < r < and r > 0. This family 
of solutions belong to the static plane- symmetric Einstein-Maxwell fields for which the 
surface element of the surfaces S2 with metric e~'^^WQAB has a non- vanishing gradient 
(see Chapter 15.4 in [7j). It can also be regarded as the fiat counterpart of the Reissner- 
Nordstrom metric. Although that family of spacetimes in [7] presents, in principle, two 
parameters m and e, whenever m 7^ a convenient change in r can bring both e and m 
in [7j into a single parameter. If m = that family falls into the Rab = case. 
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5.2.2. Case (ii) This case is characterised by a ^ab with positive constant curvature 
—4AB = 4\AB\. Coordinates {^9, can therefore be chosen such that 

nABdx^dx^ = ^^M^^ + sin^M^^), 



where G (0, vr) and (/? G [0, 27r). After the change e^^ = and renaming a = 2\J (? — 1, 
the direct substitutions lead to the hne-element 



d.s^ = -— -dr + 



{AR^ + B 



i2 



4L45 ^ ^ 



(36) 



(i?2_l)2 a2(Ai?2 + 5)4 

for an electromagnetic potential given by 

^ = e^oaR^ 
2i?2-l- 

Note that although three parameters appear in the metric, one of them can be absorved 
applying a convenient change of coordinates, and therefore only two are relevant. Now, 
this metric contains two very different subfamilies, depending on whether the gradient 
of the surface element of the (f} surfaces (see above) vanishes or not. Direct 
computation shows that the gradient vanishes if and only ii A + B = 0. 

When A + B ^ one must obtain the Reissner-Nordstrom solution. Indeed, the 
change {t, x} — )■ {r, r} given by 



l-AJ\AB\(c'^-l)Br J\AB\(c'^-l) 

e^- = V — , r = 26 ^' -t, (37) 

1 + 4J\AB\{c^ ^l)Ar A + B 



where = 1, followed by the rearranging of the constants A, B into 

f R — A 

Qc = , M = , , (38) 

A\AB\V^^ 8\AB\y^V^^ ^ ' 

leads to the Reissner-Nordstrom metric in canonical coordinates 



ds' = -\l- — + %\dT^+(l- — + %\ \r^ + r^id'd' + sin'Mip''), (39) 



in the ranges < r < M — yM^ — Q1 and M + yM^ — Q1 < and its corresponding 
electromagnetic potential 

r 

after a trivial shift. Note that — > by construction (see below) and that the 
usual Q and P [8] obviously correspond to cos 6Qc and sin 6Qc respectively. 

The line-element of the special family for which A + B = {] can be conveniently 
writen as 



ds' = - sinh^ {^)dT^ + dz" + 6' {d^^ + sin^ M^^) (40) 



for 



$ = e^" cosh ( ^ 
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after the changes r = a/c^ — It and sinhx = [sinh(?//6)]~^, where 1/6 = AA^y/c"^ — 1. 
This is the well known Bertotti-Robinson solution, which is also characterised by 
being the only homogeneous Einstein-Maxwell field with a homogeneous non-null 
Maxwell field, and the only conformally fiat solution with a non-null Maxwell field 
[7]. Furthermore, the Bertotti-Robinson solution is known to describe the near-horizon 
limit of an extreme Reissner- Nordstrom black hole [21] . 

It is worth noticing here that the relationship 4:\AB\{M'^ — Ql) = {A + B^Ql 
implies that in this class (ii) of solutions we are only finding the — > part of 
the Reissner-Nordstrom solution. Indeed, the extreme case — = is excluded 
in this class (ii) because A + B = in (l36l) leads to the Bertotti-Robinson solution 
instead. This is due to the fact that in this case (ii) we are considering solutions with 
R > whereas the extreme Reissner-Nordstrom solution has R = 0, thus falling into 
the Majumdar-Papapetrou class. The — < case implies R < 0, and therefore 
will appear in the case < 1 below. To sum up, the line-element ( 136|1 corresponds to 
the (static and — Ql > 0) Reissner-Nordstrom solution containing the near-horizon 
Bertotti-Robinson metric as a limit instead of the extreme case. Note, again, that only 
two parameters in (l36l) are relevant, but for the sake of shortness we do not pursue the 
rewritting of ( 136|) any further. 



5.3. Case (Hi) 

This case is characterised by a Qab with negative constant curvature —4AB. 
Coordinates {i}, ip} can therefore be chosen such that 

QABdx^dx^ = -^{d^^ + sinh2i9ci(^2)^ 

where {} G (0, oo) and (p G (—00,00). As in the previous case (ii), after the change 
e^^ = R^ and a = 2a/c^ — 1, the direct substitutions lead to the same line-element ( !36l) 
with sin'i? changed by sinh^?. 

Since AB > in this case, A + B cannot vanish, and therefore the change (1371) 
is always possible. After performing the same parameter redefinitions ( 1381) one obtains 
the metric 

ds' = -(-l-^ + 9l]dT'+(-l- — + %] \r' + r'idd' + sinh^ ^V)- (41) 
y r J \ ^ ^2 j 

and its corresponding electromagnetic potential 

r 

after a trivial shift. In this case the only constraint on the values of the parameters M 
and Qc is + ^ 0. The range for the coordinate r for which the metric is static is 
given by — M — ^ fvP + Ql < r < —M+ \J M'^ + Ql- This is the hyperbolic counterpart 
of the Reissner-Nordstrom solution. 
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5.4- Case < 1 

The equation that differs from the previous case > 1 is 

Rab = -2VaVb. (42) 

We proceed in an analogous way to solve the system (l42l) . (I3T]) and fl32l) . The difference 



in sign in ( H2l) compared to ( 1301) only affects the equation for W, whose solution is given 
now by 

ly = (Ae" + Ze-*^)-2, (43) 

where A is a complex number. The same previous procedure shows now that the surfaces 
(5*2, ^ab) are of positive constant curvature AAA. Coordinates (f} can therefore be 
chosen such that 

flABdx^dx^ = -^{d^^ + siii^M^^), 

where d G (0,7r) and G [0,27r). The complete line-element of the solution is found 
using ( 143|) on (133|) and taking into account (l26l) for V = x. This case is analogous to 
the case (ii) above. When A + A 7^ 0, as expected, the change of coordinates 



^2^^ ^ ^-^^AA{\-c^)AT ^ ^ ^ jAA{\-c^) ^ 
2 + 4y^AA(l-c2)Ar' ^ + ^ 

and the renaming 

i(A- A) 



^AA^fY^' %{AAfl'^^fV^ 

is what takes us to the Reissner- Nordstrom metric ( 139|) . but for — < 0. Note 
that with the above definitions AAA[M'^ - Q^) = -{A + A)'^Q'^. 

If A+A = the change z = acotx with = -(A-A)^^! - and T = Vl - cH, 
which induces the change — — c^, leads to 

ds^ = - 1^1 + dT^ +(^1 + ^^ dz^ + a^d'd^ + sin2 Mip^), (44) 

and the electromagnetic potential (after a trivial shift) $ = —e^^z/a. The metric 
corresponds again to the near- horizon Bertotti- Robinson spacetime (HOj) . now in different 
coordinates {T, z}. 

Let us stress the fact that the "intrinsic" difference that has led to fHOj) and fjUl) 
in the present setting lies in the different sign of the scalar curvature R of the scaled 
quotient space hat with respect to the Killing vectors dr and dx, respectively, but it 
is not an intrinsic property of the spacetime. In other words, the difference lies in the 
possibility of choosing timelike Killing vector fields in the Bertotti-Robinson spacetime 
with associated positive and negative curved scaled quotient spaces hat- Note, however, 
that in the Reissner-Nordstrom case the dt Killing is intrinsically defined (unit at 
infinity) and that the sign of R corresponds to the sign of — Q^, which leads to 
two globally different spacetimes. 
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Table 1. Possible quotient spaces (E3,/iab) in conformastat electrovacuum solutions. 



hab Ricci scalar 








^ = 




Majumdar-Papapctrou 




^ > 


flat 


Plane-symmetric fields 






spherical 


Bertotti- Robinson 








Reissner-Nordstrom exterior 






hyperbolic 


hyperbolic Reissner-Nordstrom 




^ < 


=>spherical 


Bertotti- Robinson 








Reissner-Nordstrom exterior 


M2 - Q2 < 



6. Results 

The combination of the above theorems and the classification of the functionally 
dependent conformastat electrovacuum solutions in Section [5] leads to the following 
final result: 

Theorem 3 All conformastat electrovacuum spacetimes either belong to the Majumdar- 
Papapetrou class or correspond to either 

the Bertotti- Robins on conformally flat solutions [JW ; 

the non- extreme exterior Reissner-Nordstrom solution (3^) . 

or its flat ( f^) or hyperbolic ( fTip counterparts. 

■ 

Let us stress that the five classes are exclusive, and that the extreme Reissner-Nordstrom 
case is included in the Majumdar-Papapetrou class. For completeness we include the 
Table [1] with a classification of the conformastat electrovacuum solutions in terms of 
the geometrical properties of the timelike static congruence defined by dt in ([T]) with 

Aa = 0. 

The first corollary of this theorem and the classification presented in Table [1] 
constitutes an improved local characterisation of Majumdar-Papapetrou. The original 
local characterisation (see e.g. [7]) states that it is the class of static electrovacuum 
spacetimes with flat hab- Here we have relaxed the requirement on hab by showing that 

Corollary 3.1 The Majumdar-Papapetrou class of solutions are the static electrovac- 
uum spacetimes with conformally flat hab ^i^d R = 0. 

An alternative statement of the above theorem is that the static charged black hole 
related geometries, that is, the Majumdar-Papapetrou, the Reissner-Nordstrom exterior 
and the near-horizon Bertotti-Roinson geometry, together with the trivial plane and 
hyperbolic generalisations of Reissner-Nordstrom, are locally characterised by being the 
only conformastat electrovacuum spacetimes. A global argument regarding asymptotic 
flatness can then be used to establish that 
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Corollary 3.2 The conformastat electrovacuum asymptotically flat spacetimes are 
either isometric 

• to the asymptotically fiat subset of the Majumdar-Papapetrou class 

• or to the exterior Reissner- Nordstrom solution. 

Further global considerations may be finally used to single out the black hole 
geometries whithin the Majumdar-Papapetrou class, the so-called standard Majumdar- 
Papapetrou, favoured by the uniqueness results in [1]. In order to do that one should ask 
for the global requirements that single out the standard Majumdar-Papapetrou among 
the complete class that appear as hypotheses in the results shown in [25], which basically 
consist of demanding a non-empty black hole region and a non-singular domain of outer 
communications. 

Acknowledgments 

We are grateful to Marc Mars for his many suggestions, criticisms and comments on 
this work. We also thank Jose M M Senovilla for his suggestions and careful reading of 
this manuscript. All the computations shown in the appendices have been carried out 
by programming the whole procedures in REDUCE. GAG wants to thank the warm 
hospitality of Dept. de Fisica Teorica, UPV/EHU, where this work was made. RV is 
funded by project IT-221-07 from the Basque Government, and thanks support from 
project F1S2007-61800 (MICINN). 

Appendix A. Proof of /^s = 

In this Appendix we present the proof of Proposition [1], as indicated in Section S) the 
proof that the integrability conditions of the equations fl2Tl) . fl22|) and fl23l) imply an 
additional isometry. This follows, exactly up to a couple of points and modulo some 
typos and missing terms in intermediate steps. Sections 4 and 5 of [15]. Let us recall 
that the two differencies of our proof with that in [15] come simply from the two aspects 
in which the treatment of the static electrovacuum case differs to that of the stationary 
vacuum case, as explained in Section [31 

The first is the fact that our functions a and 7 = 5 are two composed functions, 
one conjugate to the other, and not one complex function and its complex conjugate. 
The same goes for the coordinates = a and = a. Although the product a'-f must 
be positive (see Subection 14. ip . other factors such as a27i(= ^252) can be negative in 
general. The positiveness of Q!27i(= a2tt2) in the complex case in used precisely in the 
final step of the proof in [15], Section 5. Therefore we will need some further steps to 
complete the proof in our case. 

The second difference comes from the number N (see ([H])), which infers a different 
numeric factor in one composed equation. This difference will only imply different 
combinations to produce the equations needed in each step of the proof. We will indicate 
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all the calculations keeping an arbitrary A^. The purpose is twofold. Apart from the 
usual completeness reason, we also want to reproduce the proof in [15], and by doing so, 
indicate (and fix) some intermediate errors (typos and some missing factors) we have 
found in [15], Section 5. Therefore we will keep using the notation ~ for the conjugate 
operation that particularises to the complex conjugate in the complex case. 

The starting point is the set of equations fl2Tl) . fl22|) and the equations for the Ricci 
tensor fl23l) . Note that only enters one equation in f l23|) . the (1,2) component. The 
aim is to prove that /3 does not depend on x^. To do so, we assume P3 ^ in order 
to find a contradiction. Recall that neither a nor 7 depend on x^. The first step is 
to strictly follow the arguments in [15], Section 3, where the integrability conditions 
for the functions /3 and p in the equations (|2T|) . (l22l) are obtained. The integrability 
conditions eventually yield three differential equations, namely (29b), (29c) and (29d) 
in [15], together with their conjugates, for the functions a and 7. 

The second step follows Section 4 in [15], in which the equations for the Ricci 
tensor ( 1231) are used. The generalisation to include an arbitrary A^ is straightforward 
and we simply indicate the equation involved and the result. A^ appears in the 
R12 = N/{x^ + x^)^ equation component, and therefore contributes (only) to equation 
(31) -with (32)- in [15], which now reads 

- + P"^^-)^- + + ^R.) = 4^^;;^^. (A.1) 

This equation (and its conjugate) is convenient because, after using the equations for the 
derivatives of p and /3 (equations (24) and (25) in [15]), provides the only combination 
in which no d^a derivatives appear, leading to a polynomial of degree 9 in /3. The 10 
coefficients of the polynomial must thus vanish, providing, in principle, 10 differential 
equations for a and 7. Nevertheless, those 10 equations are proportional to two 
independent composed equations plus one imaginary equation. Indeed, a straigforward 
calculation shows that the equations corresponding to the odd powers of (3 are all 
multiples of the composed equation 

{(T + afli2 - 4(a + ahi - 2{2N + 1)7 = 0. (A.2) 

A^ only affects the odd coefficients, and thus this is in fact the only equation where A^ 
appears. The equations for the even powers of f3 provide the composed equation 

(a + af{a-fu + 57722 - 3ai7i - 373) + 127^ + {a + a)(14a7i - 672) = (A.3) 

plus the imaginary equation 

{a + a)2(722 - an) - 6((t + a){j2 - ai) - 12(a - 7) = 0. (A.4) 

Equation flA.2l) particularises to (33) in [15] for A^ = 1, and flA.3|l and flA.4l) correspond 
to (34) and (29b) in [15] respectively. As claimed in [15], the composed equation (1A.3P 
implies (29d) in [I5] and one can easily check that ( 1A.2I) implies (29c) in [15]. As stated 
in [15], there may appear another combination of the equations for the Ricci tensor, 

namely aRu — 'jRii{= aRu — 7-R22) = 0. However, this equation provides no new 
information. All in all we are finally left with equations ( 1A.2I) . ( 1A.3I) and ( ]A.4I) . 
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Appendix A.l. The system of PDEs for a and 7 

Summing up, the complete system of equations for a and 7 which decouples from the 
rest of the field equations is given by (]A.2p . (1A.3I) and ( 1A.4I) . which are conveniently 
rewritten as 



2 Ml 

^i = -7i2 + -7i + y;:^ = (A.5) 



3 3 

E2 = -722 + an + -(72 - tti) + ^(a - 7) = (A. 6) 



1 3 

E3 = a7ii + 57722 - 3(ai7i + 73) + -(7a7i - 3772) + ^7^ = 0, (A. 7) 

where r = (cr + cf)/2 and M = 2N + 1. Since we will be only interested in M = 3 and 
M = 9 we will implicitly assume at some points that certain polynomials in M with 
other roots do not vanish, and in fact, that M > 0. The 'accent is used here to keep an 
analogous notation to that in [15] , and the only purpose is to denote differently certain 
equations. Note, however, that the ' here corresponds to the tilde in [TB] . Note that 
E2 = —E2, and therefore the above system of equations contains 5 real equations. 

The procedure consists of generating new differential equations by computing the 
integrability conditions of the system {Ei, Ei, E2, E3, E3). This procedure will be fixed 
by the use of very specific sets of rules, which must be applied in strict order. Before 
setting the rules, let us produce two useful combinations after using Ei and Ei to 
eliminate 712 and au respectively: 

£■3 = 57E2 + -E3 = "711 + 57aii - 3(ai7i + 73) 

+ -(7a7i + 12772 - 157ai) + ^(15a - I27) = 0, (A.8) 
/2 \ 2 ' 1 

^4 = (92 -^3 + -7 - 672 -£2 E3 = ^2711 - 72^22 + -(27ttll + 3^172 + a27i) 

\r / r r 

1 /2M + 3 2M + 9 \ / ^«7 x 

+ 7"i - «7i - 3«72) - (2M + 9)-J = 0. (A.9) 

The first and main rule is 

(i) Multiplication by unknown functions (or their derivatives) is allowed only when the 
resulting equation does not exceed the cubic degree in the unkown functions. 

This rule only affects the choice of combinations to generate new equations. Since we are 
going to indicate these combinations explicitly, this rule does not need to be implemented 
in the algorithm. It must also be stressed that in all the equations the factors that will 
be isolated (and thence "eliminated" ) appear linearly and with a non-zero multiplicative 
factor. The first set of rules, as such, Rulesi = {{ii), {Hi), (iv), (f), (vi)} reads: 

(ii) eliminate 712 and au using Ei and Ei respectively. 
[Hi) eliminate the product a27ii by using i?4. 
{iv) eliminate the product 7^22 by using E3. 
(v) ehminate the product 0711 by using E^. 
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iyi) eliminate 722 by using E2. 

[vii) eliminate the product 71^22 by using rule {vi) applied to £"4. 

In what follows we simply indicate the chain of equations used, and the explicit 
expressions will be only given when needed. For the sake of concreteness we prefer 
to specify whenever any set of rules is applied to any expression / by Rules{f). 
The sequence of equations starts with 

= Rulesi[ 7(92 -£4 ], 

and follows with 

Eg = Rulesi[ r^d^E^ - 2rE5 - Ua^E^i ], 

3 ' 3 ~ ~ 
Efs = Rulesil 3^ _ 2M ^^ ^ 2^^^^ " "1^3) ], 

2 ' 4 ' ~ 30 ~ 

E7 = r^ g^gg_^^ i?n/egi[ 292^6 - -Ee + (M - 18)d2{Rulesi[^iE3]) - —71^3 ]• 

Note that in the third factor, as indicated, one must apply some rules before 
differentiating. The chain of equations follows with 

rlo^ 7M + 27 ^ (5M + 18)(2M-51)~ 1 , 
Es ^ Rules,[ -d2E, ^a2E, + ^ ^^^m) " r^^ 1' 

which results in a first order equation. From this point onwards it is convenient to define 
a new set of rules (keeping the first rule (i)): Rules2 = {{ii), (vii), (Hi), (iv), (vi)}. The 
chain follows with 

Eg = Rules2[ diEs ], 

and a new sequence given by 

F5 = Rules2[ 71^^2-5^4 ], 

Fe = Rules2[ ^a^s - -F5 ], 

r 

Eg = Rules2[ d2Es ], 

F-r = Rules2[d2EQ--Ee], 

r 

used to construct 

Eio = Rules2[ Ad2Eg + d2{r^FG) - 2rFe - ^Eg - ^Fg ] 

Eu = Rules2\ 2d2(r^Eio) - ArEw + 4(13M + 56) (-Eg + -E^ 

V r 

(2M2-M-141)(2M-51)~ M + 2,_. _ ,1 ^ 
^ '-F7 ]—{M^ + 10M + 5)—aE3 



6(M-33) ' 4 ' r2 



(9M3 + 106M2 + 645M + 876) 1 ^ 

\ 7TTT 7:E' 



4(M + 5) r 



2' 



^7 



The general explicit expression of Eu reads 

/ 2 \ ^ 

-24r2(M + 5)Eu = Ui - -aj 7 (322M^ + 3771M2 + 12018M + 7713) 
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+2 (ai - 7ia (319M^ + 4560M^ + 19473M + 24888) 
/ 2 \ ^ 

+ (72 - -7 j a (148M3 + 7797M2 + 64416M + 147303) 
+2 (12 - -1) "27 (346M3 + 3831 + 14262Af + 20433 



r 

The procedure follows by taking the imaginary part 

4 

3(33-My(M + l)' 

from where 



E12 = Rules2[ diEn ], 
F12 = Rules2[ d2En ], 

^ „ , , 2/0^ 2,^ ^ A 23Af2 - 236M - 483 ^ 
Ei3 = Rules2[ 2r^ f ^2^12 - 7(^12 + ^^12) ) , ^1 

ii^i3 reads, explicitly, 



r' -V 9(M 



Ll 



^ ^ 27(M + 1)^ ^ 

^ 8(55M3 + 1187M2 + 6087M + 8091) 

7 (^ai - - a ^72 - ^7^ = 0. 

The next equation is given by 

1,^ 29M + 141 ^ , / 2 \ / 2 \ 

Ea = -{En - 3(^^i) -g^) = (^72 - -7j «27 " [ai - -a) 7i« = 0. 

Let us recall here that the only positiveness property we can use in the general case 
is a7 > 0. 

Let us set up a new set of rules Rules ab = {{b), (a)}, where 
[b) eliminate the factor a7| using Ei,, 
(a) eliminate the factor 02727 using Ea, 
which applied to En leads to 

r 



r\M + ^)RuleSab[ ] = ( oi - -a 



-a + 2i?7ia 
r 



(A.IO) 



with 

B = — f-665A'/^ - 8391M2 - 33735M - 45321 
24 V 

A = (-235M^ - 5784M2 - 38217M - 77508) . (A.ll) 

Note that A < and i? < since M > 0. The factor {^ai — ^a) cannot vanish, 
since otherwise Ei would lead to a = 0, and the same argument holds for the factor 
(72 — ^7) using El and 7 7^ 0. As a result, (lA.lOP and its conjugate lead to the next 
pair of equations: 

Eq = Aa (^72 - ^7) + 25^27 = 0. 
We now use this equation to set up the next set of rules RuIcsq = {(Oi), (Oii)} where 
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(Oi) eliminate 072 using Eq, 
(on) eliminate •yai using Eq. 
The chain of equations follows with 

2 2 

Eu = Rules2\ diEo Eq 1, -Fm = Rules2 \ i92-E'o -£^0 

from where we get 

Eli = — Rulesol a'^Eii ], 
which explicitly reads 

ii4 = 2a27i5 (2A~^5 + 1) + 07^ f ^AA^ + A + 5iv) = 0. (A.12) 



V2 

Note that £'14 = i?i4. Since 07 > this equation implies 0271 < 0. What we will really 
need later is simply 0271 7^ 0. From ( ]A.12I) we set up the new rule Rulesu = {(14z)}, 
where 

(14z) eliminate 0271 using £"14. 
The next equation reads 

Fi5 = Ruleso[ Fu ], 
which explicitly reads 

Fi5 = ana{A - 105) + 6alB + a272(^ + 8B) 

+ - {3aia{2B -A)- 1a2i{A + 8fi)) + S^a^ {A - IB) = 0, 
which, since A — 105 > 0, we use to set up the last rule Rulesi^ = {(15?)}, where 
(15?) eliminate aua using F15. 

The final step consists on using the previous Ea, differentiate it, and use the sets of rules 
we have just defined in a very specific order. The precise algorithm starts with 

EiQ = •yRulesi^l o:Rulesi[ (92 -^a ] ]• 

Note that at this point we have ignored rule (i), but the outcome will be precisely the 
desired result, because 

E final = Rulesu[ Rulesoi Rulesu[ RuleSab[ -Eie ] ] ] ] 
reads, explicitly, 

3 1 4((M-22)A-25M)fi2 + (3^M + 6A + 105M + 165)A2 _ 
Ef^nai - 7ia 7^ A{A + 2B){IQB-A) ~ °- 

The last factor, after using f lA.lip to introduce the values of A and B in terms of M, is a 
fraction containing polynomials in M in which all the coefficients are positive numbers. 
Therefore, the only solution to E final = would be 710^7 = 0, which is not allowed by 
virtue of (IA.12P and 07 > 0. 

We have thus shown that 07 > ^ /^s = for any positive A^, and in particular, 
in the stationary vacuum case (A^ = 1), recovering the result in [15], and in the static 
electrovacuum case (A^ = 4). 
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The starting point is equations ( |2T1) . (1221) and ( l23l) wlien all functions depend only on 
a and a. Equation f l22|) is used to isolate the two derivatives of p, 

2D^=/3/3,i-7A2 + -(7-/3)/3, 

r 



P 

2D^ = 1313,2 - oil3,i + -{a - f3)/3, 
p r 



(B.l) 



where r = (cr + (x)/2. These two equations, which are of course related by conjugation, 
will be used to eliminate p in what follows. The integrability condition will be dealt 
with later. 

From equation ( 12T]) and taking into account that 07 > 0, we can now isolate 0:2, 
and its conjugate 71, to obtain 



"2 = 7' 



-2f3ai + 072 + af3i + (3f32 + - (3a(3 - 2a-f - f3 
-2/372 + 7^1 + 7/32 + /3/3i + - (37/3 - 2a7 - /3^) 



71 



a 



(B.2) 



We use these two expressions to compute the second derivatives 012, 022, 7ii, 722 in terms 
of «!, 72, 7i2, ttii, and the first and second derivatives of (3. We use their substitutions 
in what follows. 

We concentrate now on the Ricci equations f l23|) . From the equation R^^ = we 
isolate /3i2, which reads 



2Df3i2 = 7(«i + /32)/32 + a(72 + /3i)/3i - /3(ai/3i + 72/32 + 2/3i/32) 
+ ^[(/3 - 7)«i + (/3 - a)72] + ^(a7 + /32 - /3a - /37). 
Now, from the real combination aRu + 7-R22 = we isolate /3ii, which yields 
aD(3u = 7(a72 - /3ai)/32 - /3[a(/3i + 72)72 + (7"! - 2/372)^1] + (2^7 - /3')ai/3i 



(B.3) 



(7 - /3)a(3/372 - 27/32) + 7/3(2a/3i - 3/3ai) + {f3^ - 3a-f)af3i + {2a-f + f3^)f3ai 



a 



7(5/3' 



a7 - 2/37) --/3(a7 + 3/3') 



(B.4) 



The complex cojugate equation provides /322, which solves in turn the imaginary equation 
aRii — 7-R22 = 0. This equation is in fact equivalent to the compatibility condition of the 
above system (IB.ip for p. It is straightforward to check that the compatibility condition 
/3ii2 = /3i2i is automatically satisfied. The first important consequence of (IB. 41) is that 
if /3 = the equation reduces to 7 = 0, which contradicts 017 > 0. We must therefore 
take /3 7^ in what follows. 

The equation R13 = is identically satisfied, so it only remains to consider the 
equation R12 = N/{4r^). It is convenient first to substitute ai and 72 by two new 
functions Z and W = Z defined by the relations 

[2/3/32 + 3a/3i + -a/3 --f3^ + W 
4/5 L r r . 
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2/3/3i + 37/32 + ^7/3 - + ^ 



(B.5) 



from where we will also obtain an and 712 in terms of Zi and Wi. Let us also use the 
substitutions 



/3 



/3 



bi = /3i- -, h2 = P2- -■ 
The equation R12 = N/ (4r^) thus leads to the real relation 
F = 9a7 I a6? + 760 - 20h^h2 -^(N + 2)d \ - aZ^ - -^W^ + 2I3W Z = 0. 



9r2 



(B.6) 



Appendix B.l. Case A: 



Let us assume Z 7^ 0, so that + 7!^^ — 2(3 ZW) > (this is a positive definite 

product due to l^D = ^^(07 — (3^) > 0). The procedure now consists on finding new 
equations on Z, W, Zi Wi, bi, and a,f3,'y by first differentiating F and follow by 
using the substitutions above. In the following expressions we also use ( ]B.6p by isolating 
From the first two derivatives diF and 92-F one can isolate Zi and Wi and find the 
explicit expressions: 

24a7/3D (Wi--W 
\ r 

3a7 { [(9a7 - 12I3^)W - a/SZ] bi + [(5^7 - A(3^)Z + 3-f/3W)] 62} 

+ (9a7 - 6/3^)71^2 - (14a7 - 12^^)(3WZ + (5a7 - Ql3^)aZ^ 

aWP2{aZ - m)Z + ?>ai{Zh2-Wh,) 
+ 4^10 ivjiy ^2 «Z2 + 7iy2_2/3Ziy ' ^ 

= (4a7 + QP'^)WZ - 3-f(3W^ - 7a^Z^ 



24a(3D Zi 



3a[37/3iy + (9a7 - lQl3^)Z]bi + 3[(lla7 - Al3^)l3Z - (5a7 - 2/3^)-fW]b2 



+ 4(16 - iV)«,J^^°^^'-W-3.^Z>. 



3a7(2/3Z - 71^)62 



(B.8) 



We continue by taking the 82 derivatives of (IB. 51) and use the above expressions to 
obtain Z2 and W2 in terms of Z, W, hi, 62 and a,/3,7. One can therefore investigate 
the compatibility condition Z12 = Z21 (equivalent to W12 = W21), which provides one 
real equation: 

(16 - N){-fW - l3Z)D^KFjj = 0, 

where 

K = aZ'^ -2(3ZW + -fW'^, 

and 



'II 



607 



K + 



16 -N 



+K 



a-fl3D 
D 



^ha{(3Z - -fW)] + 



K{2f3^ - a-f) - af3j- {3{N + 8)a7 + 2{N + 2)f3^) 



-4 ((14 - N)N + 32) — aV/S^ 



(B.9) 
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For the cases we are interested in we can assume 7^ 16, and since 'yW — /3Z 7^ as 
otherwise i^' = 0, we necessarily have Fjj = 0. In the complex case (stationary vacuum) 
studied in [161 [I7j one resorts to the fact that /3 = {a + af{y -y) > and iV = 1 to 
establish that K + ^(16 — N)a(3'yD > 0. In the general case, however, one must still 
consider two subcases. 

Subcase Al: K + ^(16 — N)af3'^/D 7^ 0. We use fIB.QP to isolate hi and consider 
the imaginary combination '^[{/3Z — '-fW)adi]Fjj = 0, which leads to 

a(3-f{-fW - pZ)DK^Gi = 0, 

where the factor Gj satisfies 

\Gi - WFii = -a{-fW - I3Z) (k + ~ ^ af3-fD) 2Q{3a-fb2Z - aZ^). 

\ ry ^ J 

(Note that for this last step one must not use Fu explicitly and leave bi unsubstituted.) 
As a result, the equation 

Fjjj = 3(30762^ - = (B.IO) 

follows. We proceed with a twin combination to the previous, ^[{aZ — l3W)adi\Fii = 0, 
to finish recovering the two derivatives of Fjj. This combination, after neglecting non- 
vanishing terms, leads to Gm = 0, where Gjjj satisfies 

1 / 16 - \ 
Gin - ZFii = 6a-f{-fW - 13 Z) K + ^a/^T^* %aZbi), 



V iy'2i 



and therefore yields to 

'^{aZbi) = 0. (B.ll) 

If -fW^ - aZV equations (iBlOj) and (iRU]) lead to bi = W/{3a) (and 62 = Z/{3-f)), 
which substituted on F = (1B.6P implies Da(3'-f{N + 2) = 0, which is impossible in the 
present case. Since we are interested in the cases A^ = 1 and A^ = 4 we will also assume 
in the following that A^ + 2 7^ 0. Therefore we need 

-fW^ -aZ"^ = (B.12) 

to make (IB. 101) and (]B.11|) linearly dependent. 

We continue by taking the imaginary combination {Zd2 — Wd^) (\B.12\i and 
substituting 62 from (IB. lip . Using (IB. 120 . and after neglecting non- vanishing terms, 
that combination is shown to lead to the following equation 



+ 



(15a7 - 3^^)aZ^ + (9/3^ - 21a-f)(3WZ + 6(16 - N)a'^-f'^/3 ^ 

D 



(5a7 - 9l3'){aZ' - PW Z) + 4(16 - N)a'YP 



abi 

W = Q. (B.13) 



The real combination {Zd2 + 14^(91 ) (IB. 121) is proportional to 7!^^ — aZ"^ and therefore 
bears no information. 

On the other hand, let us take (9i(lKTn| and apply the following chain of 
substitutions: first Zi and Wi from (IB.SP and (]B.7p . followed by 02 and 71 from (lB.2p . 
then use ( IB. 51) and follow by first substituting the first derivatives of /3 by bi and 62 
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and then the first derivatives of bi and 62 by the corresponding expressions in terms 
of Z,W,bi,b2,a,'-f,/3 that come from the above relations for /3i2, (and (^^22)- Next 
substitute 62 from ( IB.IOI) and use the combination of ( IB. 61) with ( IB. Ill) so that b^ can 
be isolated in terms of Z, W, a, 7, /3 only. Finally, use (]B.12|) to eliminate first the factor 
'-fW'^ and then to express the equation in the form fbi + gW = for some factors / 
and g not depending on bi, just like equation flB.131) . At this point it is convenient to 
introduce the definition 

K = aZ'^ + -fW^ - 2I3ZW, 

so that K = 2Z{aZ — f3W) because of (]B.12p and use it within the factors / and g 
to express first Z^W^ in terms of and ZWK and then Z^ and (separately) in 
terms of K and WZ. Using this procedure the expression for 9i(IKT0]) can be cast as 



(2a7 + 

+ {(5/32-2a7)fs:2 
D 



- 2(5DZWK\ bia 
2PDZWK 



+Aa-fl3— -/< (10 - 4iV)a7 - (2 + iV)/3" -A{N + 2)l3DZW 



W = 0. 



(B.14) 



We already have the equations needed to end the proof: ( IB. 91) . ( IB. 131) and (IB. 141) . 
On top of the above defined K, we will now make use of the following extra useful 
definitions 



n 



16 - D 
—-—j^O (and^-1), S^{2 + N)ajP-, 



2 + N - ... . . 

so that in the stationary vacuum case n = 5 and in the static electrovacuum case n = 2. 
Let us stress that in the general case /3 does not have a fixed sign, and thence neither 
5 has, even for N > 2. It is only the complex case that ensures us that /3 > and 
therefore 5 < for iV > 2 (recall that aj > 0, l^D > 0). 

After using ( IB. Ill) to get rid of 62 and (IB. 121) together with the above procedure 
for expressions of the form fbi + gW = so that / and g depend on W and Z only 
through the factors K and ZW, equation (]B.9P reads 



3a'^-fK{K + n6)bi + [a-f{K + A6){K + n6) - 2/3^ K{K - 6) 
Analogously, equations (IB.lSp and (IB.14P read, respectively. 



W = 0. 



(B.15) 



{5a-f - (3'^)K - ADZWK + AnSa-f abi + {5a-f - 9(3'^)K + 8n5a-f W = 0, (B.16) 



-3 



(2a7 + I3^)K^ - 2I3DZWK] abi - {{2a-f - 5l3^)K^ + 2/30 ZWK 
-4(5 [((n - 3)a7 - Z^^) K - A^DZW] - Sn5'^a^] W = Q. 



Let us now rewrite ( IB. 61) conveniently as 



(B.17) 



(B.18) 
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Since l^K > this equation implies l^{K + 86) > 0. For + 85 = it is necessary and 
sufficient that bi = 0(= 62)- In that case, though, flB.lSP together with flB.161) lead to 
n = 4 (A^ = 8/5). Since we are not interested in that case we can assume n 7^ 4 in the 
following, so that &i 7^ and thus [."^{K + 86) > 0. 

The combination of ( IB. 151) and ( IB. 161) that cancels the terms 6bi reads 

-3abiK{DK-Al3DZW) + [-K'^D + A[2l3'^ - (n - A)a-f]6K + lQn6^a-f] W = 0.(B.19) 

The combination (lRT9l) - (lRT7l) leads to 

9ir=^a6i(a7+/3^) + [3(07 - 3(3'^)K'^ + 4(3DZW{K + 86) + 4a-fK6{n -2)]w = 0.(B.20) 

On the other hand, let us isolate WZ from ( IB. 191) . use that on ( ]B.20|) . multiply the 
result by SKabi and then use (IB.18|) to get rid of b^. Again, multiply the result by 
bi/{2W) and use (IB. 181) to get rid of bf. The resuling equation, after neglecting the 
multiplying factors W and K + 86, reads 

-3abi {K-46) [{2a-f + ^^)K + 2n6a-f] -K [2(^07 - 2/3=^)5 + (2^7 - bp'^)K\ = 0. (B.21) 

Another useful combination consists on taking ( IB. 151) . multiply it by abi/W and use 
(IRTsD to get rid of bl to get 

3abi [(a7 - 2p'^)K^ + {a^{n + 4) + 2/3^)^5 + Ana^6'^] +a^{n6+K){K+86)W = 0.(B.22) 
Now, the combination -4x (IRT51) - (Ib:2T]) +2x (Ib:221) multiplied by 1/{3K^^) leads 

to 

3abiK -{K- A6)W = 0. (B.23) 

Proceeding once more by multiplying this equation by bi and using (]B.18p to eliminate 
bl, we obtain a different relation between bi and W: 

3abi{K -A6) -{K + 86)W = 0. (B.24) 

Finally, isolating bi from the latter and substituting in (IB.23P we finally obtain 

6W\6-K) = 0, (B.25) 

which now implies 6 = K because we are assuming Z ^ and /3 7^ 0. We only need now 
to isolate bi from (1B.23P and substitute that onto ( ]B.15P using also 6 = K to obtain 

KWa-f{n + 1) = 0, 

which contradicts our assumption Z 7^ in the present case. This finishes subcase Al. 

Subcase A2: K + ^{IQ — N)al3'~fD = 0. With the above definitions this is 
K + n6 = 0. We only have to go back to equation ( IB. 91) and express it in terms of 
K to obtain 

K(3{n + 1) = 0, 

which contradicts our assumption Z 7^ in the present case. This finishes subcase A2 
and therefore case A completely. 
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Appendix B.2. Case B: 

We deal now with Z = 0. Let us take ai from ( IB.SP and use it on the first equation in 
(E2I) to get 

7^2 - + + - f 27 - ]-l3] = 0, (B.26) 



2 ' r V 2 
followed by 72 from fIB.SP to obtain 

4a2/3r- - 3/32ar + 3a/3 = 0, (B.27) 

or equivalently, (a^/3~^r^)^2 = 0. The solution is thus of the form a^l3~^r^ = ('^{x'^) for 
some analytic (or hyperbolic analytic [23]) function C{x^). From this equation we have 

which used back into (1B.26P leads first to 

The first thing this equation implies is that ( is real, and thence, by the (generalised) 
Cauchy-Riemann equations di^{() = d2Q{C)y ^23^(C) = '•^f^2^^(C)) C must be constant. 
The result 

/3 = r2(ar-2 + 6)4 

for real constants a and b thus follows. Introducing this solution together with the above 
expression for a (and 7) into flB.3P leads to ab{a + br"^) = 0, which contradicts /3 7^ 0. 
This finishes case B and therefore the proof. 
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